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Abstract 

Decay-frame Kinematics (DK) has previously been introduced as a tech- 
nique to reconstruct neutralino masses from their three-body decays to lep- 
tons. This work is an extension to the case of two-body decays through on- 
shell sleptons, with Monte Carlo simulation of LHC collisions demonstrating 
reconstruction of neutralino masses for the SPSla benchmark point. 



1 Introduction 



The Large Hadron Collider (LHC) at CERN, Geneva, is expected to provide direct 
evidence for any New Physics beyond the Standard Model (SM) at the TeV energy 
scale. The properties of these particles may shed light on the origin of electroweak 
symmetry breaking and the nature of dark matter, and give clues to a more funda- 
mental theory of Nature. To explain the deficiencies of the SM, a large variety of 
theories has been put forward, and it is only by carefully measuring the properties of 
any new particles that one will be able to discriminate between them. Chief amongst 
these properties are the masses of the new particles. 

Many of these new theories (e.g., supersymmetry and extra dimensions) contain 
a stable weakly interacting massive particle that will only be 'visible' in an LHC 
detector as missing energy. Such particles are natural dark matter candidates by 
virtue of their interactions, but pose problems for mass measurements at hadron 
colliders since one cannot in general reconstruct the kinematics of an event in which 
these particles are produced. To worsen the problem, the parton-parton interactions 
in the collider have by their very nature an unknown center of mass energy. 

In the literature a number of techniques has been developed to get around this 
problem. These fall into two general classes: those that perform a fit to or set a 
limit on masses using information from the entire event sample, and those which rely 
exclusively on events near the endpoint of a kinematic distribution. Belonging to 
the first class are "Mass-Shell Techniques" (MSTs), represented by the work done 
in [1,2], [3-5] and [6], which depend on maximizing the solvability of assumed mass- 
shell constraints in a given sample of events. This has been shown to be very effective 
if enough such constraints are available^ Also belonging to this class are techniques 
which work with extrema of a "transverse mass" variable, e.g., mx2 [8~14], and 
techniques which look at the shape of complete invariant mass distributions [15-19]. 
Much research has recently become focused in these areas. 

The second class of mass reconstruction techniques most notably includes the 
traditional kinematic endpoint method [20-28], where the endpoints of various in- 
variant mass distributions can be matched to analytical functions of the unknown 
masses, that in turn can be inverted to solve for the same masses in suitably long 
decay chains. Such methods have been studied for over a decade already, and would 
seem to have been largely explored for the simplest, and most probable, decay chains 
in popular theories such as the Minimal Supersymmetric Standard Model (MSSM). 
However, there is more to be done with events near an endpoint as demonstrated 
recently in [29]: here a Decay- frame Kinematics (DK) technique utilizes the fact 
that events at a kinematic endpoint can have exactly-known kinematics in terms of 
production angles and energies of all particles in the assumed decay chain. Events 
near an endpoint will thus have approximately known decay-frame kinematics, which 
allows one to constrain and solve for unknown masses. In [29] this was demonstrated 
for the case of neutralino three-body decays through off-shell sleptons to lepton pairs 

however, problems arise if, e.g., there are three-body decays or too many invisible decay prod- 
ucts; see [7] for further discussion. 



plus missing energy carried away by the lightest supersymmetric particle (LSP), the 
lightest neutralino, e.g. x® — > £ + £~Xv The on-shell case was deferred to a future 
work — this work. 

In the following we will demonstrate the use of DK in the case of on-shell neu- 
tralino decays, i.e. Xi ~^ — > t + t~x\-, though it should be stressed that the 
technique demonstrated can be applied to any similar cascade decay process aris- 
ing in any New Physics model. Section [2] begins with a review of the off-shell case, 
demonstrating its application at a NMSSM parameter point. Section [3] then discusses 
the main new development of this paper, the generalisation of the DK technique to 
the on-shell case, where it is found that several subtleties emerge beyond what was 
found for the comparatively simple off-shell case; we further present a Monte Carlo 
(MC) study of the mSUGRA SPSla benchmark point [30], where the DK technique 
proves quite capable of reconstructing the relevant neutralino masses from x\x\ de- 
cays. Section H] gives our conclusions. 

2 Off-Shell Decays 

We begin with a brief review of the DK technique applied to the case of neutralino 
three-body decays. For a more complete treatment, see [29]. We consider production 
of neutralino pairs in the MSSM which undergo three-body decays to electrons, 
muons, and Xi (the LSP): 

pp X - X' + *°(- e + e"x?) A*Vx?), (1) 

where X represents either a Z* or any MSSM production channel via a Higgs (H° or 
A ) or cascade from gluino/squark pair-production, while X' are SM states potentially 
produced in association, relevant in this context only for the measurement of missing 
momentum. The physical observables of interest from one event thus consist of 
four leptonic four-momenta p e ± „±, from which we may construct the usual di-lepton 
invariant masses, M ee and M w , and missing momentum in two transverse directions, 
assumed equal to the sum of the two LSPs' transverse momenta, pT vl . If we happen 
to have an event where both the invariant masses M ee and M MM are maximal, as shown 
in Fig. it will be subject to the system of constraints (hereafter we abbreviate 
rrii = m^a) 



M™ 3 * = nn-rm, (2) 

M%T = rni-m x , (3) 

PJ++PJ- = 0, (4) 

fy+Pi- = 0, (5) 

(p x +Px') T = f T (observed), (6) 



where leptonic momenta are written in the rest frame of the respective parent neu- 
tralino, and p e '± = Aip e ± and p'± = A 2 p A1 ± define the appropriate Lorentz trans- 
formations Ai ; 2 from the lab frame. This system gives ten equations for the nine 
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Figure 1: (a) Neutralino three-body decays with maximal M t 
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leptons have equal and opposite momenta while the x± (LSP) is at rest, (b) If the Xi 3 
a two-body decay via on-shell sleptons, the LSPs are no longer stationary in the neutralino 
decay frame, but have momenta K^j collinear with the leptons' momenta. 



unknowns, the velocities (3x,i and the masses m\^j, allowing us to actually over- 
constrain the masses mi^j. The j3\^ which satisfy (JH) and (jSJ), making the total 
momentum of each lepton pair zero, are uniquely given by 

fl= t £++ f ^d p 2 = . (7) 

E e + + E e - + E^- 

Now, because of the condition that M ee and are maximal, the corresponding 
Ai^, which take the e + e~x? an d systems to their respective Xij res t frames, 

also bring each Xi to rest. Thus their four-momenta in these frames must be (mi, 0), 
which, when inverse-Lorentz-transformed by A^ 2 , giving (771171,2 , ^1(^7) 1,2), have 
to agree with the observed missing momentum § T \ this matching condition along 
each transverse direction (say x and y) then gives two independent determinations 
of mi. 

m'i = , a ; ** lQ ; and m'[ = - . . . (8) 

(Pxl)i + (^7)2 (A,7)i + (^7)2 

Since we are assuming that both M ee and are precisely maximal — the 
perfect event of Fig. [T^, — we should get m\ = m" = m\ from such an event. In 
practice of course, we can only expect to find an event within some neighborhood e 
of the endpoints, M ee)W , = M™^„ ± e, in which case one finds that m' x and m" are 
offset by 0(^j2em\) from m\ [29]. One might then expect that applying ([7]) and (jHJ) 
to a sample of events near the endpoint should give a distribution of m' x and m" 
peaked near m\ with a spread determined by sample purity. 

Here, to lend further support to the generality of the above technique, let us 
demonstrate its application to the rather challenging example of an NMSSM (Next- 
to-Minimal Supersymmetric Standard Model) scenario described in [19]. This has a 
supersymmetric particle spectrum containing five neutralinos, the lightest of which is 
99% singlino, and a generic feature of the parameter space that gives the correct dark 
matter density is a significant degeneracy between the singlino and second lightest 
neutralino mass. This gives rise to copious production of soft lepton pairs with small 



invariant masses, M« < 10 GeV, from the three-body decay x[> ~~ * £ + ^~Xi- F° r more 
details on this scenario see [19]. 

We study the benchmark "Point A" of that paper, a point which has M^ ax = 
9.7 GeV and mi = 105.4 GeV, using the same MC setup and fast detector simulation 
as in [19]. For details of the simulation see also Section [3] of the present paper. To 
isolate signal events of the type (jTJ), we place the following cuts on our events: 

• Require missing transverse energy > 100 GeV. 

• Require at least two hard jets with p T > 150, 100 GeV. 

• Require four isolated leptons with flavor structure e + e~[i + fi~ and px > 7 (4) 
for e (/i). All such leptons must pass the lepton efficiency cuts employed in [19] 
modeled on full simulation results given in [31]. 

From the surviving events we construct a wedgebox plot of the di-electron versus the 
di-muon invariant mass, for a number of events equivalent to 30 fb _1 of statistics 
at the LHC The result is seen in Fig. [2^, showing a clear box-like structure at 
M eetflfJi ~ 10 GeV, the endpoint of the di-lepton invariant mass distribution for the 
X° 2 decay. 

Choosing a sampling region in a rather generous neighborhood of the endpoint, 
M u = 10 ± 4 GeV gives 0(100) events. We now apply Eqs. © and flS} to each of 
these, demanding that m\ and m'[ agree to within 20%. Although only about 20 
events survive this criterion, the resulting mi distribution can be seen in Fig. [2b to 
peak quite prominently slightly below the nominal value of mi = 105.4 GeV. The 
systematic error of the method is seen to be comparable to the estimate made earlier. 

With higher statistics, this allows us to determine the absolute LSP mass to rather 
good precision. The results for 300 fb _1 of data are shown in Fig. [2b. Here we narrow 
down our sampling region to e = 1 GeV from the wedgebox edge. With a Gaussian 
distribution we obtain a best fit value of mi = 98.2 ± 3.9 GeV with x 2 /ndf = 0.41. 
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Figure 2: (a) Wedgebox plot for 30 fb^ 1 of integrated LHC luminosity for "Point A" 
of [19]. (b) The Xi mass distribution for events sampled from the boxed region in (a). We 
see a peak in rough agreement with the nominal value of mi = 105.4 GeV. (c) Same as (b) 
but for 300 fb' 1 ande = l GeV. 



3 On-Shell Decays 



3.1 Kinematics 



Let us now continue to the main focus of this paper, i.e. the added complications 
that arise when the neutralinos decay through on-shell intermediate sleptons: 



X - X' + #X?(- e^/i^ - e + e>VxiXl)- (9) 

are maximal, as illustrated in Fig.QJ), two-body kinematics gives 
the following system of constraints: 
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(10) 

(11) 
(12) 
(13) 

(14) 

(15) 
(16) 

where we have assumed a common slepton mass m s = mg = m^. See the Appendix 
for details of the derivation. The antiparallel conditions ( fT2l) and (fl3~|) force (3\$ to be 
in the planes of the respective leptons, so there are really only four boost parameters 
to find; adding the unknown masses mx^j^ to this gives eight unknowns which can 
thus be solved for by the eight constraints f|T0l) - f|T6l) . In principle, if one were handed 
an event of the type in Fig. [lb, one could numerically apply (fID~|) - (fTBj) . scanning over 
the eight- dimensional space of unknowns for a solution. Needless to say, this is not 
the most practical approach, nor particularly enlightening as to the nature of any 
solution which might be found. 

Instead, we will proceed temporarily as if we already knew the individual p x T x , as 
opposed to just their sum (fT6|) . Consider, then, just one of the neutralino decays, say 
Xi ~ * e ± e =F — > e + e^Xv I n the x° rest frame, the leptons' three-momenta (= P[ 2 ) are 
back-to-back and collinear with the LSP's three-momentum (= f 1 ); thus, in the lab 
frame where xl has a velocity (3i, the boosted momenta Pi, P%, and j* lie in the same 
plane (see Fig. [3^)- Looking at this the other way round, (3i is the Lorentz boost 
which makes the observed lepton momenta antiparallel and fixes the magnitude of 
their sum to be Ki — which we a priori don't know at this point — i.e. constraints 
(1T21) and ([111). Necessarily, p 1 must be in the observed leptons' plane. If we choose a 
basis in this plane (p», p±) parallel/perpendiculaiE] to the total leptonic momentum 



2 A potential ambiguity in defining p± is resolved by defining it such that P\ ■ p± is positive: 

Pi - (Pi ■ P\\)p\\ 
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Figure 3: (a) When lepton pairs from the decay chain x° — > £ £ T — > ( + i~Xi have maximal 
invariant mass, their momenta P\p in the lab frame are coplanar with the LSP momentum 
P. (b) The LSP's longitudinal momentum Pl can be found by the condition that J 3 is in 
the plane of the leptons. 



P = Pi + P 2 , then the boost /3i = ((3\\, (3±) must satisfy three sets of constraints: 
1. The transformed leptonic momenta must be antiparallel: 

P[-P^ = -\P[\\P^ (17) 
where the transformed four-vectors are given in terms of the boost by 
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with /3 2 = [3 2 + /3| and 7 = (1 - (3 2 )- 1/2 . 
2. The transformed total leptonic momentum must equal K { = (Kf, K^): 
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where the components (If], if^) are also unknown at this point. 
3. The inverse-Lorentz-boosted LSP four-momentum must satisfy 
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(18) 
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(20) 



After some algebra, see the Appendix, these constraints are found to uniquely deter- 
mine the unknown boost (/3m, /3j_) in terms of the known lab frame leptonic momenta 
and unknown LSP momenta: 

^ = E(l + ax) ^ ^ = a/3 »' (21) 

where 

a^^- and s S ^ ~ ^ ~ ^ ■ 
+ p 2E{ExPl + E 2 Pl) 

Thus, knowing the leptonic momenta and missing momentum from the LSP deter- 
mines x and a, hence (3\ and all the kinematic information in the event. At first 
glance this may seem useless since we can only have knowledge of the transverse 
component of If in the lab coordinate system, fo, and there are two LSPs that 
contribute to the measured total missing momentum. 

However, given the transverse component of the LSP momentum we can in fact 
reconstruct the longitudinal component ]fr by the following trick: since Jf must lie 
in the plane of the leptons while fa is by definition in the transverse x-y plane, l^r 

— # — * 

must be of the precise size along z to bring fr + fr into the leptons' plane (see 
Fig. [3b), i.e. (fi T + fi L ) • (A x P 2 ) = 0, giving 

g _ f > Tx{P\yP2z — P\zPly) + Jf > Ty{P\zP2x ~ P\xPlz) * ^\ 



P\xPly ~ P\yPlx 

With both fir and fir known we may immediately project ^ into the basis (/V -fjjjfl 
compute a and insert into (I2TI) to solve for the boosts. We can then use Eq. (|T9|) to 
solve for k} and 

which are also related by = a(Kf — P). We see that in the limit a —>■ both 
K}' ± —>■ 0, and that also Eq. ( 12TT) correctly reduces to the off-shell result (3 = P/E. 
Finally, the LSP mass can then be found from the energy component of Eq. ( }20l) : 



mi 



(7-l)£f\ K} (j-1)/3 ± KJ- 



nw^{ 1+ ^)t^^r^ ] - hr - Kr - (24) 



which again reduces to the off-shell result of Eq. (jSj) when 0. The heavier 

neutralino mass follows from the energy component of ffT§|) and energy conservation 
in its decay, i.e. 



m % = JKf + m 2 1 + 1 E-P in P, (25) 



3 If the leptons happen to be parallel fi, remains undetermined. We ignore events with this 
pathological arrangement. 



while the slepton mass is related to mi and m, by (1101) . 

Finally, let us return to deal with the realistic situation where we know only the 
sum of the LSP momenta {p x +Px') T ■ F rom the discussion above, every assignment of 
Ft = P X T will yield a set of masses {m' 1 , m'^ m' s } which satisfy f fTUj) . ffl~2]) . and ffl4|) . 
Then the other LSP has its p x T fixed as p x T = f T — p x T , giving another set of masses 
{m'l, m", m"} which satisfy f fTTl) . f[T3"j) . and ffTBT) . Under the simplifying assumption 
that the event contains the process in with i = j, we should clearly insist that 
at least {W 1; m^} ~ {m", m"} within some error (we reserve the possibility that 
m' s ^ m's)- This, in principle, provides two constraints on our choice of the two 
components of p x T , i.e. the system ( flQ|) - (fl6l) is solved]^ 

Thus, in the end, we still have to resort to a numerical search for a solution 
to (TTUT ) (TTBT) . but this is only over the two-dimensional space of one of the LSPs' 
transverse momenta, (p x , p x ). Nevertheless, it is not at all obvious that there won't 
be multiple solutions with different {mi, m^} within a given level of tolerance — and 
when one adds to this the same caveat as in the three-body case of picking events 
within some e of the endpoint, as well as the effects of detector smearing and issues 
with backgrounds, it will have to fall to a MC simulation to test the practicality of 
the method. 

3.2 Monte Carlo Test 

We perform a Monte Carlo study by generating SUSY signal events for the SPSla 
benchmark point [30] using PYTHIA 6.413 [32], and SM background events with 
HERWIG 6 . 510 [33,34], interfaced to ALPGEN 2 . 13 [35] for the production of high jet 
multiplicities matched to parton showers and JIMMY 4.31 [36] for multiple interac- 
tions. The benchmark point is chosen mainly for the sake of comparison with results 
obtained with other mass reconstruction techniques, which we will comment more on 
in Section |H The generated events are then put through a fast simulation of a generic 
LHC detector, AcerDET-1 .0 [37], widely used to simulate analyses of high-p^ physics 
at the LHC. This incorporates such detector effects as the deposition of energy in 
calorimeter cells, and the smearing of electron, photon, muon and hadronic cluster 
energies with parameterized resolutions. The AcerDET-1 .0 isolation requirement for 
leptons is less than 10 GeV energy in a R = 0.2 cone around the lepton and a min- 
imum distance of AR = 0.4 from calorimetric clusters. The MC setup is essentially 
the same as in [19] and we refer the reader to that paper for more details. However, 
we point out that we use pt dependent lepton efficiencies based on full simulation 
studies published in [31]. 

All SUSY processes and relevant SM backgrounds are generated with a number 
of events corresponding to an integrated LHC luminosity of 300 fb _1 . The dominant 
type of neutralino pairs produced at SPSla are X2X2] we wm therefore concentrate on 
decays of the form yQ(— ¥ e+e ~Xi) xlK - > Xi)i an d results in the previous section 
can be simplified somewhat by setting i = j = 2, and in particular K = Ki = Kj. 



4 Notice that something quite interesting has happened here in that we have gotten around the 
usual four-fold ambiguity in designating 'near' and 'far' leptons in the decay chains. 



For this analysis we use the same set of cuts as for the NMSSM case in the 
previous Section, giving a signal size of roughly 470 events. Note the requirement of 
four isolated leptons with flavor structure e + e~/i + /z~ reduces most SM backgrounds 
to a negligible level [38]. Support for this in the context of a full simulation of the 
ATLAS detector is found in Higgs searches for the channel h — > ZZ* — > A£, e.g. 
discussed in [39,40]. The remaining backgrounds of any importance are found to 
be ti, Zbb and irreducible Z^*'Z. We have simulated a large sample of it events 
with up to two additional hard jets, and find no surviving events with the additional 
missing energy and jet cuts. The Zbb and Z^Z backgrounds are also expected to 
be very small after these cuts. However, because the Z mass is sufficiently far from 
the dilepton edge, any remaining events from these backgrounds do not significantly 
influence the region of interest in the di-electron versus di-muon invarant mass plane, 
shown in the wedgebox plot of Fig. 

The position of the edge of the box-like structure at M« « 75 GeV is visually 
apparent in Fig. H^, and, as shown in several studies, can be brought into precise (sub- 
GeV) agreement with the nominal value M^ ax = 77.07 GeV, by the standard study 
of the flavor-subtracted di-lepton mass distribution shown in Fig. [5b. One advantage 
of the DK technique is that we do not strictly need such precise determination of the 
edge — GeV-level will do to determine our sampling region — but we will assume 
that the edge has been measured to 76.7 ± 0.1 GeV as quoted in [24]. 

Events in a broad neighborhood of the corner of the box, M ee wt = 65 ±e GeV (see 
further comments below on sampling regions), are passed to the on-shell DK analysis 
described in the previous Section. In detail, the procedure used is the following: 

1. An event is selected if the two invariant masses M ee and M w both lie within 
the e-defined region of the wedgebox plot. 

2. A point in (fax, fr y )-space is chosen by a uniform scan of —500 GeV < fa x ,y < 
500 GeV, in 0.2 GeV steps; the point is assumed to equal the transverse momen- 
tum of the LSP accompanying the e + and e~ whose four- momenta are (E, P)i,2, 
respectively. 

3. The longitudinal component of the LSP's momentum is found from (j22]l . 

4. Components of the LSP momentum in the basis parallel/perpendicular to the 
total leptonic momentum P = Pi + P 2 are determined and used to compute 
a = P L /(P\\ + P). 

5. The boost parameters (3\\ and (3±_ are now computed from fl2T]) . and K 1 - from 
d23D. 

6. The masses m' 2 } are computed from ff24l) and ff25|) . 

5 However, for some SUSY parameter points one might have a dilepton edge very close to the Z 
mass. Though the Z background events would thus be unavoidably mixed in with signal events, 
they would in general not have solutions in the numerical procedure described in the following. This 
property of DK gives it a certain resilience in the face of backgrounds. 
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Figure 4: (a) Wedgebox plot for 300 flr x of integrated LHC luminosity at the SPSla 
benchmark point. Events are sampled from the boxed region shown for DK analysis, (b) 
The flavor- subtracted dilepton invariant mass distribution provides a clean determination 
of the edge MJ^ ax = 77.1 GeV. The plot is shown for an integrated luminosity of 30 fb^ 1 . 



7. Using the missing momentum constraint p^ =f T — p x T ', steps 3-6 are repeated 
for the LSP accompanying the fi + fi~ pair, obtaining either another set of masses 
{m", m'2} or no valid solution for the second mass set. 

8. If no valid second set of masses was obtained, the point is assigned zero weight. 
Otherwise, the two sets of mass solutions are plotted with the following weight: 

P(P P )- 1 exof K-Q 2 K-^H m) 

where a is our expectation for the spread between the mass values on either side 
of the event. This should be of the order of the missing energy resolution, which 
is a function of the total transverse energy Et deposited in the calorimeters. 
We therefore assign sigma on an event-by-event basis, using the expected per- 
formance of the ATLAS detector in SUSY events (see figure 10.84 of reference 
[31]): 



a = 0.57^^2 Et. (27) 

9. The scan is continued until all points have been assigned a weight. 

10. The procedure is repeated for all events that have invariant masses sufficiently 
close to the endpoint. 



Our final mass distribution for a single event is obtained by histogramming all 
mass solutions found in the scan over missing momentum components, weighted by 




Figure 5: Mass distribution for 300 fb^ 1 of data at SPSla, obtained using the procedure 
described in the text. The contribution of individual events with degenerate solutions is 
still clearly visible, although it is a region close to the nominal mass values of (7711,7712) = 
(96.1,176.8) GeV that emerges with the largest total weight. 

Eq. (|26|) . For events that lie exactly at the endpoint, and where the mismeasurement 
of lepton momenta and missing energy is negligible, the peak of the two mass dis- 
tributions for {m'i, m' 2 } and {m'(, m 2 '} should coincide, agreeing with the nominal 
values of the masses. We find that this remains accurate for events that lie sufficiently 
close to the endpoint, and thus with a value of e that is not too large. 

In fact, choosing a value of e is essentially a trade-off between accumulating 
statistics by allowing more events to pass the cut, thus reducing the fluctuations 
that come from the smearing of lepton and missing energy momenta, and protecting 
the integrity of the mass solutions that are obtained by the procedure at e = 0. We 
find that a compromise value for e of 15 GeV gives just enough events with peak 
region close to the nominal values to dominate over events that display either one 
or two degenerate peaks or peaks in the wrong place. In the latter case, we observe 
that the maximum weight at the peak is lower. 

By summing the distributions for events within our choice of e, the resulting total 
mass distribution for 300 fb _1 of data is shown in Fig. [5J Although the existence 
of multiple peaks is clear with the limited statistics, the point with the largest total 
weight is found to be (7/11,7712) = (114.75,191.75) GeV, close to the nominal value 
of (mi,m 2 ) = (96.1,176.8) GeV. We emphasize that this procedure is simply a 
suggestion for an estimator of the masses, and although there are similarities in 



shape, Eq. (j2"E]) is not a likelihood function. 

To check the robustness of the estimator and find the statistical error on making 
such a measurement, we have performed 10 independent 'experiments' with 300 fb _1 
integrated luminosity. In each case, the mass solution with the largest total weight 
fell near the nominal masses, with a standard deviation of 20.2 GeV on mi and 
21.2 GeV on m 2 . While these errors are quite large, there is undoubtedly scope for 
improvement. By better understanding the properties of events with degenerate or 
wrong solutions one could search for a system of kinematic cuts to remove these 
subsets; one could also increase statistics to tighten the cut on e or investigate other 
estimators for the masses with better properties with respect to these events. 

Incidentally, had we wrongly assumed off-shell kinematics at this parameter point, 
and hence tried using the technique of Section [2] to analyze events in the boxed 
region of Fig. H^, we would have failed since essentially no events provide two mass 
solutions with a near equal mass. This provides a way of distinguishing a sample of 
on- versus off-shell decays^] distinct from the usual way of measuring the departure 
of the di-lepton mass distribution from a triangular shape [19], or looking for specific 
relationships between the positions of lepton-jet invariant mass maxima [41]. 

Finally, note that this method, although we have not explicitly demonstrated it, 
is in principle also applicable to extracting the slepton mass. 

4 Conclusions 

This paper completes the demonstration of the application of the DK technique to 
neutralino decays in SUSY models, or indeed any similar decay chain in models such 
as e.g. UED [42], having explicitly shown the procedure for reconstructing both off- 
and on-shell decays to lepton pairs in realistic MC simulations of two very different 
supersymmetry benchmark points. We find that in the three-body decay scenario we 
can reconstruct the LSP mass with an accuracy of around 4 GeV, while the decay 
through an on-shell slepton allows a precision of 20 GeV. 

DK has the advantage of simplicity and robustness in the face of backgrounds; for 
non-signal events there tends to be no solution for the constraint equations, or at least 
no preferred solution in the distribution of events passing loose constraint require- 
ments. Moreover, since the method only makes use of unlike di-leptons (e + e~/i + /i~) 
from the \2 ~^ Xi transitions, it is insensitive to the combinatoric issues that arise 
when one considers particles produced further up the decay chain. It may hence 
complement mass determinations which exploit the full chains such as [5,6]. Gen- 
eralization to other decaying states (e.g. charginos as in [43]), perhaps using jets 
instead of or in addition to leptons, is open to investigation. 

Perhaps the major disadvantage of DK is the requirement of a high event rate: 
one typically needs at least O(10) events in the neighborhood of a kinematic endpoint 

6 One exception occurs if it happens that K = 0, i.e. m s — ^111^2. Then both on-shell and 
off-shell techniques are equally applicable and should reconstruct the same LSP mass. At SPSla, 
K w 18 GeV, which is fairly small compared to the maximum K one can get with the same neutralino 
masses, K max = ^(1 — mf/m|) ~ 62 GeV, but safely away from zero. 



to make the reconstruction stable, and, in the case of four-lepton final states, that 
translates to several hundreds of events needed on a wedgebox plot. In addition, the 
case of on-shell decays gives rise to extra solutions in the mass space that can be 
eliminated with more statistics, but ultimately contribute to an increased error in 
the reconstructed masses. 

Thus, in the case of neutralino-pair production considered in this work, DK may, 
depending on the parameter point Nature has chosen, perhaps only serve useful as a 
check on results obtained with other techniques that do not depend on events near an 
endpoint. As mentioned in the Introduction, these most prominently include MSTs 
and rriT2 techniques, which at the SPSla point happen to work quite well. 
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Appendix 



On- Shell Kinematics 

Consider the neutralino decay 



Xi 



(28) 



in the neutralino rest frame. If the di-lepton invariant mass is maximal, all the decay 
products must be collinear (say along x). In particular, four-momentum conservation 
forces 
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where m s is the slepton mass and the lepton is assumed to be massless. Similarly, in 
the slepton's decay frame, 
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which, when boosted back to the Xi res t frame using (3 = —{mf — m 2 s )/{rril + m|), 
becomes 




/ 



V 



rrii 



2m 



-mi 



\ 

2m'i / 



and 



V 



m^—m^m^ 
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From these equations it is easy to verify (fit)]) . (1TT1) . f|T4|) . and (TT3T) . 



Finding the Lorentz Boost Parameters 

Here we indicate in more detail how one may arrive at Eq. (1211) . Starting from 
Eq. (1201) we have two relevant equations, 

n+Kl = P\\ (ly/x? + "»? - (7 - 1)^J - (7 - !)^) , 
+ ^ = #l (7^ + - (7 - 1)^J - (7 - 1)^) • 
Taking the ratio of these one obtains after some rearranging 



(29) 



This can now be inserted into two of the equations of (fT51) . 



Kf = -/? ||T £ + P + P( 7 -1)JL, (30) 



K = -P±JE + (j - 1)^P, (31) 



which can then be used to solve for the ratio 

° = W = PT^' (32) 

Starting from the antiparallel condition (TT7T) and expanding with the lepton momenta 
from ( Tl8|) . one arrives, after some algebra, at 

= (1 + p\ + /3i)£i£ 2 - 2^||EiP| - 2p x E 1 P 2 ± ~ W\\E2P\ - 2f3 ± E 2 P^ + 
2/5||/3 ± (Pj l P 2 ± + pjp^) + (1 + /3f - PI)P}P} + (1 - P\ + Pl)P 1 ± P 2 ± , 

using that E\ j2 = | -Pi,2 1 - Substituting for a and /5y from (1321 . we get a quadratic 
equation for 

= a/3j - 26/3|| + c, 

where 

a = P!P 2 + pj'pj - P^P^ + 2a(Pj l P 2 ± + pJp^) + a 2 (PiP 2 - Pj'pJ + P^Pf), 
6 = PiPj + E 2 p} + a(P 1 P 2 ± + P 2 P 1 ± ), 

c = PiP 2 + pj'pj + P^P£. 

Using again that Ef 2 = pj' 2 H-Pf 1 ^ 2 , one can show that b 2 — Aac = 0, so this quadratic 
has a double root 

g = E l Pl + P 2 P| + ajE.Pj- + E 2 Pf) 

" PiP 2 + P|P| - Pj-Pj- + 2a(PfP 2 ± + P^Pj) + a\E x E 2 - Pj'pj + Pj-P£-) ' 

Making repeated use of that, by definition, P^ = —P 2 and P = P," +p|, one arrives 
at the equation for /3\\ in ( |2TI) . 
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